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ABSTRACT 

We study the dimensionless spin parameter j (= cJ/(GM 2 )) of uniformly rotating neutron stars and 
quark stars in general relativity. We show numerically that the maximum value of the spin parameter 
of a neutron star rotating at the Keplerian frequency is j max ^0.7 for a wide class of realistic equations 
of state. This upper bound is insensitive to the mass of the neutron star if the mass of the star is 
larger than about 1 Mq. On the other hand, the spin parameter of a quark star modeled by the MIT 
bag model can be larger than unity and does not have a universal upper bound. Its value also depends 
strongly on the bag constant and the mass of the star. Astrophysical implications of our finding will 
be discussed. 

Subject headings: Dense matter — stars: neutron — stars: rotation 



1. INTRODUCTION 

The general stationary vacuum solution (the Kerr 
spacetime) of the Einstein equations is specified uniquely 
by the gravitation al ma s s M a nd the angular momen- 
tum J (see, e.g., IWaldl ((l98l L If J < GM 2 /c, we 
have a rotating black hole. However, if J > GM 2 /c, the 
Kerr spacetime would have a naked singularity without a 
horizon. One could then consider closed timelike curve s 
and causality would be violated (|Chandrasekharl IT983) . 
While its validity has not yet been proven, the cosmic- 
censorship conjecture (|Penroselll969l) asserts that naked 
singularities cannot be formed via the gravitational col- 
lapse of a body. For this reason, it is believed that astro- 
physical black holes should satisfy the Kerr bound j < 1, 
where j = cJ/GM 2 is the dimensionless spin parameter. 

While the value of the spin parameter j plays a funda- 
mental role in black-hole physics, it appears that this is 
not the case for other stellar objects. In particular, there 
is no theoretical constraint on the value of j for stars. 
It is known that the spin parameter of main-sequence 
stars depends sensitively on the stellar mas s and can be 
much large r than unity (|Kraftlll968l 119701: iDickd 11970b 
lGravl ll~982L On the other hand, the spin parameter of 
compact stars has not been studied in detail (see below). 
As we shall discuss in more detail in Section [3J the spin 
parameter of compact stars is interesting in its own right 
for two reasons. (1) It plays a role in our understand- 
ing of the observed quasi-periodic oscillations (QPOs) in 
disk-accreting compact-star systems. (2) It determines 
the final fate of the collapse of a rot ating compact star. 

Ever since the seminal work of HartH (| 19671 ) who 
considered the limit of slow rotation, rotating compact 
stars have been studied extensively in general relativ- 
ity. In the past two decades, various different numer- 
ical codes have been developed to construct rapidly 
rotating stellar models in general relativity. We re- 
fer the reader to IStergiouIaa ((2003) for a review. As 
the rotation frequency / is a directly measurable quan- 
tity for pulsars, it is thus reasonable that the max- 
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imum value for / (i.e., the Keplerian frequency /k) 
has been one of the most studied physi cal quanti- 
ties fo r relativistic rot at ing st ar s (see, e.g.. ICook et al.l 
(fl99l: IHaensel et al l (fl995h: iKoranda et al.1 (| 19971 )7 
iBenhar et al.1 (|2005( ): IHaensel et al.1 (|2009[ )). However, 
in contrast to these previous works, we shall focus exten- 
sively on the spin parameter j. 

It has been known that the spin parameter for a 
maximum-mass neutron star lies in the range j ~ 
0.6 — 0.7 fo r most realis tic equations o f state (EOS; e.g., 
ICook et all (|1994[ ) and ISalgado et all (|1994h ). In this 
work, we first extend the previous works and show that 
the upper bound for the spin parameter j max ~ 0.7 is 
essentially independent of the mass of the neutron star if 
the gravitational mass of the star is larger than ~ 1 Mq . 
Next we study the spin parameter of self-b ound quark 
stars, which was not consid ered p reviously in lCook et al.l 
(|1994| ) and ISalgado et al.l (|1994[ ). We find that the be- 
havior of the spin parameters of neutron stars and quark 
stars is very different. In contrast to the case of neutron 
stars, the spin parameter of quark stars does not have 
a universal upper bound. It also depends sensitively on 
the parameter of the quark matter EOS and the mass of 
the star. Furthermore, the spin parameter of quark stars 
can be larger than unity. This leads us to propose that 
the spin parameter could be a useful indicator to identify 
rapidly rotating quark stars. 

The plan of this paper is as follows. Section [2] presents 
the main numerical results of this work. In Section G2 we 
discuss the astrophysical implications of our results. Our 
conclusions are summarized in Section [4] 

2. NUMERICAL RESULTS 

2.1. Numerical method and EOS models 

We make use of the numerical code rot star from the 
C++ LORENE library 1 to calculate uniformly rotat- 
ing compact star models in general rel ativity. The code 
uses a multi-domain spectral method (jBonazzol a et all 
1998) to solve the Einstein equations in a station- 
ary and axisymmetric spacetime with a matter source 

1 http://www.lorene.obspm.fr/ 
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Fig. 1. — Gravitational mass as a function of the central energy 
density for non-rotating compact stars constructed with the chosen 
EOS models in this work. 
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Fig. 2. — Keplerian frequency as a function of the gravitational 
mass for rotating compact stars constructed with the chosen EOS 
models in this work. 



(|Bonazzola et al.l Il993t iGourgoulhon et al.lll999t ). The 
code has been tested extensi vely and compared w ith a 
few different numerical codes (|Nozawa "eTatlfl99cl . 

As theoretical calculations for dense matter at 
supranuclear densities are poorly constrained, the EOS 
in the high-density c ore of compa c t star s is not well un- 
dersto od (see, e.g., IWeber et al.l (|2007|) ; lHaensel et al.l 
lj2007h for reviews). In this work, we employ eight 
realistic nuclear ma tter EOS to model rotating neu- 
tron stars: model A (|Pandharipanddll97lD . model APR 
(lAkmal et al.|[l998h. m odel AU (the AV14+UVII model 
i n IWiringa et al.1 (|1988f ) is joined to iNeeele fc Vautherinl 
1973h). model BBB2 (fBaldo et all I1997D. model FP S 



Pandharipande fc Ravenhalll 119891: lLorenz et al.lll993t ). 



model SLY4 (jDouchin fc Hacnscl 2000), model UU (the 
UV14+UVII model in IWiringa et~aTl (|1988t ) is joined 
to INeeele fc Vautherinl (Tl973l)) and mod el WS (the 
UV14+TNI model in IWiringa et all (|1988h is joined to 
lLorenz et al.1 ()1993f )). For the quark star models, we use 
the simplest MIT bag model wit h non-interacting mass- 
less quarks (jChodos et a l. 1974). Two different values, 
60 MeV fm~ 3 and 90 MeV fm" 3 , are chosen for the bag 
constant B. These values correspond approximately to 
the range of B within which the h ypothesis of strange 
matter is valid ([Haensel et al.ll2007|) . 

To illustrate the diversity of the EOS models used in 
this work, we plot the gravitational mass M against cen- 
tral energy density p c for non-rotating stars constructed 
with the chosen EOS models in Figure [TJ The quark star 
models QMB60 and QMB90 in the figure correspond to 
the cases B = 60 MeV fm~ 3 and 90 MeV far 3 , respec- 
tively. The maximum mass of compact stars depends 
quite sensitively on the EOS models and it ranges from 
about 1.5 M© to 2.2 M Q . In Figure EH we plot the Kep- 
lerian frequency /k against the gravitational mass M of 
rotating compact stars based on our chosen EOS models. 
Similar to the maximum mass for non-rotating compact 
stars, Figure [5] shows clearly that /k depends strongly 
on the EOS model. It is also sensitive to the mass of 
the star. This is the well-known reason why searching 
for rapidly rotating compact stars can provide us con- 
straints on the EOS of dense matter. 
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Fig. 3. — Maximum spin parameter as a function of the gravi- 
tational mass for rotating neutron stars. 



2.2. Neutron stars 

Now we turn to the main focus of this work: the dimen- 
sionless spin parameter j. Having seen that /k depends 
strongly on the EOS and the mass of the star, it may be 
quite surprising to learn that the maximum value of the 
spin parameter j max (as set by the Kepler limit) is quite 
universal for rotating neutron stars. In Figure [31 we plot 
jmax against the gravitational mass M for the selected 
nuclear matter EOS models. In the figure, each line cor- 
responds to one particular EOS and each point on a line 
corresponds to a star model with a fixed M rotating at 
its Keplerian frequency. Note that each sequence in the 
figure is terminated at the stellar model with the same 
total particle number as the stable maximum-mass non- 
rotating configuration. In the figure, we see that j ma x 
lies in a narrow range ~ 0.65 — 0.7 for the eight different 
nuclear matter EOS models. In particular, the values 
do not depend sensitivel y on the mass of the star. This 
exten ds previous works (|Cook et al.lll994t ISalgado et al.l 
11 9941 ) which focus on maximum-mass neutron star mod- 
els. 

While the spin parameter of an astrophysical black hole 
is constrained by j < 1, we find that the spin parameter 
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Fig. 4. — Spin parameter is plotted against the scaled rotation 
frequency for neutron stars constructed with the FPS EOS. Each 
line represents a sequence of fixed total particle number. Each 
sequence is labeled by the gravitational mass of its non-rotating 
configuration. 

of a neutron star is bounded by j max ~ 0.7. The upper 
bound jmax is quite universal for different EOS models 
and gravitational mass larger than ~ 1 Mq. For lower 
mass neutron stars, M < 1 Mq, we find that j decreases 
with decreasing M. However, we shall only focus on 
mass M > 1 Mq in this work as the observed masses of 
neutron stars are typically larger than 1 Mq. We refer 
the reader to ISteiner et al.l ()2010f ) for a recent review on 
the observed masses of neutron stars. 

So far we have studied the maximum spin parameter 
imax of neutron stars rotating at their Keplerian frequen- 
cies /k- However, realistic neutron stars in general ro- 
tate slower with frequencies / < /k- Is the spin pa- 
rameter j still insensitive to the EOS and mass of the 
star? In Figure 21 we plot the spin parameter j against 
the scaled rotation frequency ///k for the FPS EOS. In 
the figure, each line represents a sequence of fixed total 
particle number (the so-called evolutionary sequence). 
Each sequence is labeled by the gravitational mass of 
its non-rotating configuration Mq. The solid line corre- 
sponds to Mq = 0.94 Mq, the dashed line corresponds to 
Mq = 1.43 Mq, and the dashed-dotted line corresponds 
to Mq = 1.73 Mq. Note, however, that the gravita- 
tional mass increases with rotation frequency. The figure 
shows that the spin parameter changes only by at most 
10% (depending on ///k) when M changes from 0.94 
to 1.73 Mq. It is also interesting to note that the differ- 
ences between the curves decrease as the scaled frequency 
///k tends to 1 . This is quite different from the case of 
quark stars which will be studied shortly. In Figure El 
we plot j against ///k for three different EOS models. 
The stellar models on the three sequences have the same 
total particle numbers such that their baryonic masses 
are fixed at Mb = 1.6 Mq. It is seen that, for a fixed 
scaled frequency, the spin parameter of neutron stars is 
essentially independent of the EOS models. 

2.3. Quark stars 

Let us now consider self-bound quark stars using the 
MIT bag model. In Figure [6l we plot j max against M for 
the two quark matter models QMB60 and QMB90. In 
contrast to the case of neutron stars (see Figure [3]), we 
see that j ma x depends sensitively on the mass of the star. 
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Fig. 5. — Spin parameter is plotted against the scaled rotation 
frequency for neutron stars constructed with three different EOS 
models. The baryonic mass is fixed at Mb = 1.6 Mq. 

For the QMB60 model, j max is decreased by about 24% 
as M increases from 1 Mq to 2 Mq. Comparing to the 
case of neutron stars, it is also seen that j max has a more 
significant dependence on the EOS parameter, namely 
the bag constant in the MIT bag model. Analogous to 
Figure [4] for neutron stars, we plot j against ///k for 
the QMB60 model in Figure [7] Each line represents a 
sequence of fixed total particle number. As in Figure HI 
each sequence is labeled by the gravitational mass of its 
non-rotating configuration Mq. We show three different 
sequences in the figure: Mq = 0.82 Mq (solid), 1.27 Mq 
(dashed), and 1.55 Mq (dash-dotted). We see that the 
differences between the curves increase significantly as 
the scaled frequency increases. At ///k = 1, the spin 
parameter is increased by about 27% when Mq changes 
from 1.55 to 0.82 Mq. 

We see that the spin parameter of quark stars can be 
significantly larger than the upper bound (j max ~ 0.7) 
for neutron stars. It can also be larger than the Kerr 
bound j = 1 for black holes. This suggests that the 
spin parameter could be a useful indicator to identify 
rapidly rotating quark stars. Discovering even one single 
compact star with spin parameter j > 0.7 will provide 
a strong evidence for the existence of quark stars. Fi- 
nally, it should be pointed out that the fact that the 
spin parameter of quark stars can be larger than 0.7 i s 
also evident from the results of lStergioulas et ail (|1999h . 
Using the values of the gravitational mass M and angu- 
lar momen tum J for Kepleri a n qua rk stars presented in 
Table 1 of IStergioulas et ail (|1999D . it is easy to check 
that the quark stars considered by the authors all have 
jmax > 0.7. In particular, their results also suggest that 
jmax decreases with increasing M as we have seen in Fig- 
ure [6] 

3. ASTROPHYSICAL IMPLICATIONS 

We have studied the spin parameter of uniformly ro- 
tating compact stars in general relativity. Our numer- 
ical results show that the behavior of the spin parame- 
ter of quark stars is quite different from that of neutron 
stars. In particular, the spin parameter of neutron stars 
is bounded above by j max ~ 0.7, while quark stars can 
have a value larger than unity. In this section, we shall 
discuss (in our view) the astrophysical implications of 
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Fig. 6. — Maximum spin parameter as a function of the gravi- 
tational mass for rotating quark stars. 



our results. 

First, how could the spin parameter of a compact star 
be measured? Unfortunately, so far there is no gen- 
eral technique to infer the spin parameter j of compact 
stars directly. As far as we are aware, the spin param- 
eter of a compact star could be potentially measured in 
disk-accreting compact-star systems. In particular, the 
neutron stars (or quark stars) in low-mass X-ray bina- 
ries (LMXBs) provide the most natural cosmic labora- 
tories for studying the spin parameter. In order to un- 
derstand how the spin parameter might be inferred in 
disk-accreting systems, it should be noted that the spin 
parameter of the central compact star directly affects the 
particle motion around the star. For example, to first 
order in j, the orbital frequency of a point particle in a 
prog rade orbit around a compact star is given by (see, 
e.g., Ivan der Klisl (|2006t l) 



^ \ rc 2 J 



3/2 



(GM\ 



1/2 



(1) 



where r is the orbital radius. For infinitesimally tilted 
and eccentric orbits, the disk particles will have radial 
(v r ) and vertic al (i/g) epicycl i c freq uencies which also de- 
pend on j (see Ivan der Klisl (|2006f ) for the expressions). 
Furthermore, the combination vg — v r also gives rise to 
the periastron frequency of the orbit. These frequencies 
in general depend on M and j, and hence their obser- 
vations (possibly needed to be combined with the mea- 
surement of other stellar parameters such as the mass) 
would in principle provide useful information on the spin 
parameter. In fact, there are strong evidences that these 
frequencies have already been observed in LMXBs. 

It should be noted that existing algebraic relations, 
such as Equation (TTJ), which relate various orbital fre- 
quencies to stellar parameters are in general only valid 
for small spin rate j ~ 0.1. The difficulty of obtaining 
the corresponding algebraic relations for rapidly rotat- 
ing stars (j ~ 1) lies in the fact that there is no exact 
analytic representation of the vacuum spacetime outside 
a rapidly rotating compact star. Having such an ana- 
lytic representation of the spacetime metric will allow 
one to obtain the desired algebraic relations for a rapidly 
rotating compact star. A starting point along this di- 
rection would be to take the closed-form asymptotically 




Fig. 7. — Spin parameter is plotted against the scaled rotation 
frequency for quark stars constructed with the QMB60 EOS. Each 
line represents a sequence of fixed total particle number. Each 
sequence is labeled by the gravitational mass of its non-rotating 
configuration. 



fla t solution of t he E instein-Maxwell system obtained 
by iManko et al.l (|2000t ) and study the geodesies in this 
spacetime. If there is no charge and magnetic moment, 
this analytic solution depends only on the mass, angu- 
lar momentum, and quadrupole moment of the space- 
time. Furthermore, this solution only involves rational 
functions, which helps to simplify the ana l ytical study 
of geodesic motions. iBerti fc Stergioulasl J2004) have 
demonstrated that this analytic solution can describe 
the exterior spacetime of a rapidly rotating neutron star 
(e.g., j > 0.5) very well. Nevertheless, further investiga- 
tion is needed to check whether this analytic solution is 
also valid for rapidly rotating quark stars. 

One of the well-observed features of LMXBs is the 
high-frequency (~ kHz) QPOs. To date, QPOs have 
been observed in more than 20 LMXBs. These QPOs 
often come in pairs with frequencies v u and v\. In all 
systems in which the spin frequencies of the compact 
stars JAjtar have been measured, the frequency separa- 
tion Av = v u — vi is approxi mately equal to v s t a.r or 
^star/2. We refer the read er to Ivan der Klisl (12006ft and 
lLamb fc Boutloukosl ((2008) for recent reviews. While the 
physical mechanism responsible for producing the high- 
frequency QPOs is not known yet, most physical models 
involve orbital motion and disk oscillations. Hence, the 
frequencies is r , u$, v§ and various combinations of them 
are often invoked (either directly or indirectly) to explain 
high-frequency QPOs. This is the reason why one might 
hope to obtain useful information on the spin parameter 
of the central compact star in an LMXB by observing 
its QPOs. For example, if the higher frequency of the 
QPO pair (v u ) is identified with the orbital frequency 
(V(t>), one can then obtain a relation between the spin 
parameter and an upper bou nd on the mass of the com- 
pact star (Mill er et al.| [l998'). On the other hand, based 
on the so-called relativisti c precession model of QPOs 
UStella fc VietrilfliM [l999h , iTorok et al.1 (I20T0I) have re- 
cently derived a constraint relating the mass and spin 
parameter of the compact star in Cir X-l. These works 
thus demonstrate the possibility of measuring the value 
of j with an independent determination of the mass and 
vice versa. 



5 



The above brief review of high-frequency QPOs serves 
to point out the astrophysical relevance of the spin pa- 
rameter and in what situations could it be potentially 
measured. Now we are ready to discuss how one could 
make use of the finding in this work to obtain useful in- 
formation regarding the central compact star in LMXBs. 
Suppose a single well established model for the high- 
frequency QPOs can be agreed upon in the future, then 
it is possible that rapidly rotating quark stars could be 
identified from the inferred spin parameters. As shown 
in our numerical results, if the inferred spin parameter of 
the central star is larger than ~ 0.7, then the star could 
be a quark star. Of course, the star can either be a neu- 
tron star or quark star if the spin parameter is less than 
0.7. 

On the other hand, for the present situation where 
many models are available, our finding could still be 
used to put constraints on the physical models for QPOs. 
For example, let us consider the resona ntly excited disk- 
oscillation model for QPOs proposed bv lKatol (|2008l ). In 
the model, Kato suggests that the high-frequency QPOs 
are inertial-acoustic oscillations on a deformed disk that 
are resonantly excited by nonlinear couplings between 
the oscillation modes and the disk deformation. Kato 
applies the model to Cir X-l and finds that it describes 
the observed QPOs quite well if the mass of the cen- 
tral star is about 1.5 — 2.0 Mq and the spin parameter 
is j ~ 0.8. However, our work suggests that uniformly 
rotating neutron stars cannot have j > 0.7. If Kato's 
model is correct, then our finding implies that the cen- 
tral star in Cir X-l could be a quark star. On the other 
hand, if other measurements in the future (such as the 
mass-radius relation or cooling property of the central 
star) suggest that the compact star in the system is in- 
deed a traditional neutron star, then our finding would 
rule out the resonantly excited disk-oscillation model. 

After discussing how one might use the spin parameter 
to distinguish between neutron stars and quark stars in 
LMXBs, let us now turn to a different issue concern- 
ing the collapse of a rotating star to black hole. In 
the past decade, general relativistic simulations of ro- 
tational collapse of neutron stars modeled by the poly- 
tropic EOS ha ve been performed (Shibata ct al.l 12000a ; 
lShibata|[2003aHbl : iDuez et al.H2004tlBaiotti et al.ll2005aD . 
These simulations show that no stable massive disks can 
be formed around the resulting black holes. The implica- 
tion is that the collapse of a uniformly rotating neutron 
star (with j < 1) to black hole cannot lead to the black- 
hole accretion model of gamma-ray burst, which requires 
a rota t ing bl ack hole surrounded by an accretion disk (see 
iPiranl (|2005l ) for a review). It is expected that the initial 
spin parameter of the collapsing star must be j > 1 in 
order to form a massive disk aroun d the final black hole 
(|Shibatall2Q03aUbl : IDuez et al.ll2004D . As rapidly rotating 
quark stars can have j > 1, it is thus possible that the 
collapse of a rapidly rotating quark star could be a pro- 
genitor for the black-hole accretion model of gamma-ray 
burst. Furthermore, it is noted that the collapse of mas- 
sive stars to black holes (the so-called collapsars) could 
form massive disks around the black holes and hence pro- 
duce long gamma-ray bursts. Since quark stars have 
smaller mass (~ M©) comparing to massive stars, the 
collapse of quark stars might only produce small disks. 
As the duration of the accretion, and hence the timescale 



of the bursts, depends on the mass of the disk, the col- 
lapse of quark stars might thus be a mechanism for short 
gamma-ray bursts. 

On the other hand, what if the collapse process does 
not lead to any (or little) mass ejection as in the col - 
lapse of neutron stars? In fact, iBauswein et al.1 (2009) 
have recently simulated quark star mergers based on the 
MIT bag model and the conformally flat approximation 
to general relativity. Their results show that the mass 
ejection depends sensitively on the bag constant B. In 
particular, they find that there are some binary mergers 
without mass ejection. Could this conclusion, namely 
the absence of mass ejection for some values of B, also be 
true for the rotational collapse of quark stars? For the 
collapse of rotating neutron stars to black holes, apart 
from the small amount of total mass energy M and an- 
gular momentum J carried away by gravitational radi- 
ation (B aiotti et al.l l2005bl ). the final black holes have 
essentially the same M and J, and hence the same spin 
parameter j, as the initial star. For a rapidly rotating 
quark star with initial spin parameter j > 1, if there is 
no mass ejection, how could the spin parameter be re- 
duced efficiently in order to form a regular black hole 
that satisfies the Kerr bound j < 1 at the end of the 
collapse? If a Kerr black hole could not be formed in the 
process, then what would be the final fate of the collapse? 
These questions deserve further investigation using fully 
general relativistic modeling. The hope is that studying 
the collapse of quark stars might lead to the discovery of 
some new phenomena which are not seen in the collapse 
of neutron stars. This might then help to shed more 
light on the highly nonlinear dynamics of gravitational 
collapse in general relativity. 

The astrophysical implications discussed above depend 
crucially on the existence of rapidly rotating quark stars 
with j > 0.7. However, it is known that rapidly ro- 
tating compact stars may be subject to different kinds 
of secul ar or dynamical non-axisymmetric mode insta- 
bilities (Andersson 120031 ). One might thus worry that 
rotating quark stars (if they exist) may be limited 
to small spin rates due to various instabilities, and 
hence rendering a possible distinction between neutron 
stars and quark stars unlike ly to happen. However, 
IGondek -Rosih ska et al.l (|2003| ) have shown that, taking 
into account realistic values of shear viscosity, viscosity- 
driven bar mode instability cannot develop in quark stars 
modeled by the MIT bag model in any astrophysically 
relevant temperature windows. Even taking the unre- 
alistic assumption of infinite shear viscosity, the insta- 
bility can develop only if the ratio of the rotational ki- 
netic energy to the absolute value of the gravitational 
potential energy T/W is larger than 0.1375. The ex- 
act value depends on the stellar mass. For compar- 
ison, a quark star modeled by the MIT bag model, 
with B — 60 MeV fixT 3 , M = 1.146 M , and j = 
0.8, has the value T/W = 0.1 26. We have followed 
the sa me numerical procedure of IGon dek-Rosihska et all 
(2003) to check that this quark star is indeed sta- 
ble against t he viscosity-driven instability. W e refer 
the reader to [ Gondck-Rosihska fc Gourgoulhonl (2002) 
and IGondek-Rosihska et al.l ((2003) for more details 
on the numerical pro c edure . On the other hand, 
IGonde k-Rosih ska et al.l (|2003l ) have also discussed that 
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the gravitational-radiation driven r-mode instability 
seems to be unimportant for quark stars in LMXBs if the 
strange quark mass is m s c 2 ~ 200 MeV (standard value) 
or higher. The r-mode instability can develop and may 
limit quark stars to small spin rates only if the strange 
quark mass takes the relatively low value m s c 2 ~ 100 
MeV. 

Even if the viscosity-driven and gravitational- 
radiation-driven instabilities (which are both secular ef- 
fects) cannot develop in a quark star, the star may still 
subject to dynamical bar-mode instability which occurs 
at a higher spin rate. This kind of instability has not 
been studied for rotating quark stars. However, we can 
still obtain some insight from the Newtonian theory of a 
rotating incompressible star, which is a good approxima- 
tion to quark stars because of their rather uniform den- 
sity profile. The onset of dynamical inst ability occurs at 
T/W w 0.27 for an incompressible star (|Chandrasekharl 
fl96flh . For comparison, a quark star modeled by the MIT 
bag model, with B = 60 MeV fin"" 3 , M = 1.4 Mq and 
j = 1.11, has the value T/W = 0.23. Although general 
relativity may change the critical value T/W for the on- 
set of the instability to a somewhat smaller value than 
that is suggested by the Newtonian theory 2 , it is still 
very likely that rapidly rotating quark stars with j > 0.7 
can exist. Some of them may even break the Kerr bound 
j = 1 for black holes. 

4. CONCLUSIONS 

In this paper, we have studied the dimensionless spin 
parameter j of uniformly rotating neutron stars and 
quark stars in general relativity. We find that the max- 
imum value of the spin parameter (as set by the Kepler 
limit) of neutron stars is bounded above by j max ~ 0.7. 
This upper bound is essentially independent of the EOS 
of the neutron star. It is also insensitive to the mass 



of the star if the mass of the star is larger than about 
1 Mq. On the other hand, the spin parameter of quark 
stars behaves quite differently. We find that the spin pa- 
rameter of quark stars modeled by the MIT bag model 
can be larger than unity. It also depends sensitively on 
the EOS parameter (i.e., the bag constant) and the mass 
of the star. 

We have discussed (in our view) the astrophysical im- 
plications of our finding in detail in Section |3l We have 
discussed how the spin parameter of compact stars could 
be potentially measured in LMXBs and its relevance in 
the physical models for high-frequency QPOs. As a first 
application, our finding implies that the compact star in 
Cir X-l could be a quark star if the reson antly excite d 
disk-oscillation model for QPOs is correct (|Katol F2008). 
since the model requires that the spin parameter of the 
central star to be j ~ 0.8 in order to fit the observed 
QPOs. We have also speculated on how the collapse of a 
rotating quark star might be different from the collapse 
of a neutron star. As explained in Section[3l the collapse 
of a rapidly rotating quark star with j ' > 1 might lead to 
the formation of a stable disk around the resulting Kerr 
black hole. This implies that the collapse process might 
form the central engine of a gamma-ray burst. However, 
a fully general relativistic dynamical calculation (which 
has not been done in this work) is required to shed light 
on the issue. 

In conclusion, our work suggests that discovering even 
one single compact star with spin parameter j > 0.7 
will provide a strong evidence for the existence of quark 
stars, and hence verifying the hypothes is that strang e 
quark matter could be absolutely stable (|Wittenlll984 ). 
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2 Shibata ct al. (2000b) have shown that the effects of general 
relativity only change the critical value to T/W ~ 0.24 — 0.25 for 



stars modeled by a polytropic EOS. 
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